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Abstract 

We use chiral perturbation theory to investigate hadronic properties in strong electric and mag- 
netic fields. A strong-field power counting is employed, and results for pions and nucleons are 
obtained using Schwinger's proper-time method. In the limit of weak fields, we accordingly recover 
the well-known one-loop chiral perturbation theory results for the electric and magnetic polariz- 
abilities of pions and nucleons. In strong constant fields, we extend the Gell-Mann-Oakes-Renner 
relation. For the case of electric fields, we find that non-perturbative effects result in hadron decay. 
For sufficiently strong magnetic fields, the chiral analysis confirms that the nucleon hierarchy be- 
comes inverted giving rise to proton beta-decay. Properties of asymptotic expansions are explored 
by considering weak field limits. In the regime where the perturbative expansion breaks down, the 
first-order term gives the best agreement with the non-perturbative result. 
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I. INTRODUCTION 



The study of quantum field theories in the presence of external fields is a rich subject. 
In the pioneering work of Euler and Heisenberg, the one-loop effective action for QED 
was exactly computed for the special case of a constant electromagnetic field. This result 
was later formalized by Schwinger [l| in an essentially modern way. Such non-perturbative 
solutions provide important insight into vacuum structure, renormalization, and the behavior 
of pertubative approximations, see jl] for a review of effective actions and more recent 
developments. 

Probing a complicated theory like QCD with external electromagnetic fields provides a 
lever-arm with which to study the modification of vacuum and hadron structure. This in 
turn gives insight into the underlying QCD dynamics. In external electromagnetic fields, 
model field theories, and low-energy models of QCD have been studied 0, 0, 0, @, 0, @, 0] , 



as well as the QCD vacuum [10|, [TJJ. A striking fact about the QCD vacuum is that it 
spontaneously breaks chiral symmetry, which is evidenced by the non-zero value of the 
quark condensate. Because electromagnetic fields couple directly to quarks, the behavior 
of the condensate, and hence the ground state of QCD, will be affected by external fields. 
There have been numerous studies of the effects of strong fields on the condensate. Focusing 
on the model-independent investigations using chiral perturbation theory, the calculation of 



the condensate in a constant magnetic field at zero quark mass has been carried out [12l . [13 



and additionally calculations have included the effect of finite temperature [14| . The external 
field dependence of the chiral condensate has also been recently investigated in the more 



realistic scenario of non- vanishing quark mass in constant electromagnetic fields [15l. Il6l] . As 



the condensate lowers the vacuum energy, it can be expected that the addition of a magnetic 
field (which contributes positively to the action density, cx B 2 ) will strengthen spontaneous 
chiral symmetry breaking by increasing the condensate. On the other hand, the addition 
of an electric field (which contributes negatively to the action density, oc — E 2 , and has 
associated instabilities) weakens chiral symmetry breaking. These heuristic expectations 
are verified by the explicit model-independent calculations using chiral perturbation theory. 

In this work, we consider the modification of hadron structure in strong fields. We use 
chiral perturbation theory to investigate pion and nucleon properties in constant electric and 
magnetic fields. Chiral dynamics affords us the systematic and model-independent tool to 
investigate the modification of hadron structure in external fields. In electric and magnetic 
fields, we determine the effective mass of charged and neutral pions, as well as of the proton 
and neutron. The results for charged and neutral pions are shown to satisfy a generalized 
form of the Gell-Mann-Oakes-Renner relation. While the derived results are at a fixed order 
in the chiral expansion, they are non-perturbative in the strength of the external field. For 
sufficiently weak fields, the derived mass shifts of the pions and nucleons correctly repro- 
duce the well-known electric and magnetic polarizabilities from one-loop chiral perturbation 
theory. Expanding about the weak field limit allows us to investigate the behavior of per- 
turbative approximations, and directly confront the breakdown of perturbation theory. We 
see explicitly the behavior of perturbative approximations beyond the expansion's regime 
of validity. In this situation, the leading term gives the best approximation, which can 
be remarkably good. We speculate about the presence of this phenomenon in a few other 
systems. 

Our study of hadrons in strong electric and magnetic fields is additionally motivated 
by lattice gauge theory simulations [l7| • The structure of the QCD vacuum and spectrum 
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can be studied in the presence of electromagnetic fields using external field techniques with 
lattice QCD. For small field strengths, the response of hadrons is governed by their multipole 
moments and multipole polarizabilities. These are coefficients of the first few terms in a 
weak-field expansion of a hadron's two-point function in the external field. The introduction 
of constant gau ge fields on a torus, however, requires quantization conditions on the strength 



of the field 18, [191 l2(J, 21, 22, 23 . To achieve the smallest values of constant fields, the 



generic quantization condition has the form 

where a is the lattice spacing, F^ v represents the non-vanishing components of the field 
strength, the number of the lattice sites in the direction specified by fi, and n is an integer. 
Satisfying the quantization condition on current lattices, Nj = 24, inevitably leads to non- 
perturbative effects from pion loops because the ratio eF^/m^ is not considerably smaller 
than unity. Furthermore as the lattice pion masses are brought down to the physical point 
(at fixed lattice sizes), this ratio will exceed unity. Extraction of moments and polarizabilities 
will then be obscured by non-perturbative effects. To address these effects, we determine 
the behavior of pion and nucleon two-point functions in external fields to all orders in the 
ratio of field strength to pion mass squared. Furthermore the study of hadronic two-point 
functions in strong external fields extends the method proposed in [24j for extracting charged 
particle properties from lattice QCD. 

An outline of our presentation and key results are as follows. In Section [TTJ we study 
the effects of strong electric and magnetic fields on pions. To begin, we review the chiral 
Lagrangian in Section III and derive the charged pion propagator in a strong-field power 
counting. One-loop computations are pursued in Sections III Bl and III C\ where we derive 
the shift in the effective mass due to external magnetic and electric fields, respectively. 
Expanding the non-perturbative results in powers of the external field allows us to explore 
the behavior of asymptotic expansions. We extend and demonstrate low-energy theorems for 
pions in external fields in Section III Dl The nucleon is the subject of Section IHII We review 
the chiral dynamics of nucleons and deltas in Section IIII Al Effects of strong magnetic fields 
on the nucleon are investigated in Section [III Bl where spin-dependent and spin-independent 
contributions are computed at one-loop order. We confirm that the proton /3-decays in strong 
magnetic fields. Surprisingly the leading and next-to-leading contributions to this effect 
arise from terms that are completely local in the effective theory. The study of the nucleon 
in an external electric field is carried out in Section IIII CI Technical details concerning 
the computation of loop diagrams with mesons and nucleons using the strong-field power 
counting are contained in Appendix |A] and [BJ while analytic continuation is addressed in 
Appendix O In Section IIV| we end with concluding remarks, and speculations about the 
breakdown of perturbative expansions. 



II. PION IN STRONG FIELDS 
A. Chiral Lagrangian 

To investigate the properties of hadrons in strong external fields, we turn to chiral per- 
turbation theory (25[. Chiral dynamics is based solely on the symmetries of low-energy 
QCD; and, as such, provides a model-independent description of hadrons in terms of a few 
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parameters determined from phenomenology. When electromagnetically gauged, the chiral 
theory can be used to derive electromagnetic properties of hadrons. 

In the absence of external fields and at zero quark mass, the Lagrangian of QCD has a 
U(2) L <g> U(2) R chiral symmetry 1 that is reduced to SU(2) L <g> SU(2) R <g> U(1) B by the axial 
anomaly. The U(1)b symmetry corresponds to global vector symmetry for each flavor, and 
leads to baryon number conservation. The vacuum state of QCD spontaneously breaks the 
two-flavor chiral symmetry: SU(2)l®SU(2) r — > SU(2)v- The resulting Goldstone manifold 
can be parametrized by E which is an element of the coset SU(2)l <g> SU(2) R /SU(2)v, and 
the Goldstone bosons, which are the pions, are non-linearly realized in E: 

E = exp(2#//), = . (2) 

The dimensionful parameter / arises from QCD dynamics, / ~ Aqcd; and, when weak 
interactions are considered, / can be identified with the chiral limit value of the pion decay 
constant. In our conventions, we have / = 130 MeV. 

Under a chiral transformation: L £ SU(2)l, R £ SU(2)r, we have the transformation 
E — > LTjR^ . We can build an effective theory of low-energy pions by forming the most general 
chirally invariant Lagrangian. Additionally we must take into account sources of explicit 
chiral symmetry breaking. For our considerations, these are classical electromagnetic fields 
and the quark masses. The electromagnetic fields explicitly break the chiral symmetry 
down to the product of diagonal subgroups: SU(2) L <8> SU(2) R — ► U(1)l ® U(1) R . Thus 
in the absence of quark masses, only the neutral pion is a Goldstone boson. The quark 
masses explicitly break chiral symmetry down to the vector subgroup. With non-vanishing 
quark masses and non-zero electromagnetic fields, the resulting low-energy theory has only 
a U(l)y ® U(1)b symmetry. Taking into account the pattern of spontaneous and explicit 
symmetry breaking, the chiral Lagrangian has the form 

C = £ < D M EtD M E > ~ < m Q (Et + E) >, (3) 

where we have used a bracketed notation to denote flavor traces: < A >= tr(A). The quark 
mass matrix is given by itlq = diag(m u , m^), and we choose to work in the limit of strong 
isospin symmetry, m u = rrid = m. The action of the covariant derivative is specified by 

D A1 E = 9 M E + ie^[Q,E], (4) 

where Q is the quark electric charge matrix, Q = diag(2/3, —1/2), e the unit of charge, with 
e > 0, and A^ is the gauge potential for the electromagnetic field. The parameter A arises 
from QCD dynamics, A ~ Aq CD . Physically A is the value of the chiral condensate in the 
chiral limit (technically it is the negative of the condensate). 

In writing the Lagrangian, we have only included the lowest-order terms. These terms 
scale as 0(e 2 ), where e is a small dimensionless number. We assume the power counting 

■ e . (5) 



A l A l A l 



Here we consider only two light quark flavors. There is empirical evidence from lattice QCD simulations 



that three-flavor chiral expansions are ill-fated 



2G, 



27, 



28, 



23 ] . This suggests that the strange quark mass 



is possibly too large to be considered a perturbation about the SU(3) chiral limit. Stran ge h adrons then 
must be considered as heavy external flavor sources in the relevant SU(2) multiplets [30Ll3lil32|. 
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FIG. 1: Pion propagator in strong external fields. Depicted are the C(e 2 ) Born couplings which 
must be summed to derive the propagator in our power-counting scheme. 



Here A x is the chiral symmetry breaking scale, A x ~ 2^/2nf, k represents a typical pion 
momentum, is the pion mass, and F^ v is the electromagnetic field strength arising from 
the gauge potential, F^ v = d^A v — d v A^. We work with constant electromagnetic fields 
throughout. The leading-order pion mass can be determined by expanding the Lagrangian 
to quadratic order in zero field, 

f 2 ml = 4Am. (6) 

This is the Gell-Mann-Oakes-Renner relation. We will return to its validity in the presence 
of external fields below in Section III Dl Notice the counting of the field strength allows for 
relativistic effects because eF^/m 2 ~ e°. Consequently the leading-order propagator for 
charged pions is not that for a relativistic free particle. The Born couplings to the charge 
are non-perturbative. 

We will focus on two cases for the external field. The first is a magnetic field B = Bz, 
which for concreteness we take to be specified by the gauge potential 

A^ = (—5x2,0,0,0). (7) 

Of course we could always make a gauge transformation, and arrive at the same magnetic 
field. When possible, we will express our answers in terms of gauge invariant quantities. 
The second case is an electric field E = Ez. To work with the electric field, we always start 
with a Euclidean space gauge potential, for concreteness take 

A^ = (0,0,-^X4,0), (8) 

and address the analytic continuation £ — > %E. In Euclidean space, the difference be- 
tween magnetic and electric fields arises because of our choice of time direction. We will 
see this when we calculate nucleon properties using the coordinate-space LSZ reduction in 
Sections IIIIBI and IIII CI In Minkowski space, by contrast, there is considerable difference 
between magnetic and electric fields. 

As mentioned above, the neutral pion propagator maintains a simple Klein-Gordon form. 
The charged pion propagator requires summation of charge couplings, see Figure [TJ This 
summation can be done by directly considering the charged pion action Eq. ([3]) at 0(e 2 ): 

£ = D^-D^ + + ml 7i-7i + , (9) 

with = + ieQA^, where Q — ±1. Specializing to the case of the magnetic gauge 
potential Eq. ([7]), we have 

C = 7i~ (-d 2 + ml + x 2 2 e 2 B 2 + 2ix 2 eBd x ) vr + . (10) 

Using the momentum mode expansion 

r rlk 

= J ^ e "X^), (ii) 
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where k refers to the three non- vanishing components of = (ki, 0, ks, k^), the pion action 
becomes 

S = J dx 2 J ^3 7T- (x 2 ) (--^ + k 2 + ml + x\e 2 B 2 - 2x 2 eBk^j irf(x 2 ). (12) 



It will be easier to change variables: 



Between states of good k, we need to find the propagator D, which as a concrete operator is 

D- 1 = 2 (ip& + \e 2 B 2 X 2 + \[k 2 -k\ + m 2 ]) = 2(H + El/2). (14) 

The X and px are canonically conjugate. The transverse energy is defined by Ej_ = k 2 + 
/c| + m 2 . The propagator .D(X',X) is then 

1 /"°° 

L>(X',X) = (X'\D\X) = - ds (X'\e- sn \X)e- sE ^ /2 , (15) 

2 Jo 

which can be thought of as a proper-time evolution of a harmonic oscillator [l[. Given the 
form of the harmonic oscillator propagator 



(X', 0) = J CB - exp | ^— - [(X' 2 + X 2 ) cosh e5s - 2X'Xl 1 , (16) 

x ' 1 ' 1 V 2tt sinh eBs [2 sinh eBs LV ; J J ' v ; 



we have thus derived the charged pion propagator 

1 f°° , f dk ju.m^/ j h 



2J J (2vr) 3 \ z eB 



D{x',x) = - ds T^e ik < x '- x Ux' 2 -^-,s x 2 - ) e -^ /2 . (17) 



eB 



We do not include an electromagnetic gauge link. 2 The propagator in Eq. (TlTI) is a Green's 
function satisfying the relation 

(-D^D^ + ml)D(x,0) = 5\x). (18) 

This relation will be useful in simplifying calculations. 

The charged pion propagator in the case of a constant electric field, specified by Eq. flB]), 
can be written down from the magnetic result by inspection. We find 



2 Appending the Abelian gauge link, L(x, 0) = d£ M exp[ze^4 M (C)], leads to the modified Green's function 
relation 

{-d^ + ml) L(x, 0)D(x, 0) = 5 4 (x). 
In this form, we see the action of the partial derivative is equivalent to the gauge covariant derivative. 
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FIG. 2: Diagrams for a pion in an external field. Shown are the one-loop diagrams for the pion 
two-point function in an external field. The dashed lines represent the pion propagation in the 
background field. The wiggly lines terminating in crosses represent couplings of the vertices to the 
background field. 



with E\ = k\ + h\ + m^. At this level, the only difference between propagators is the 
particular labeling of coordinates. 

Based on the power counting in Eq. (151) , we need to include higher-order terms so that we 
can calculate amplitudes to one-loop level. Two further terms that depend on the external 
field strength must be added to the Lagrangian. 3 These higher-order terms scale generically 
as 0[e A \ and will only contribute at tree-level. They have the form 

C = -iotaF^ < QDpED v Tt + QD^DJL > -a 10 F^F^ < QSQS 1 " > . (20) 

These terms lead to energy shifts of charged particles proportional to the field strength 
squared. Terms with higher powers of the external field are suppressed by the chiral sym- 
metry breaking scale, and first occur in the 0(e 6 ) Lagrangian. There is an additional field- 
dependent term that can potentially contribute from the fourth-order Lagrangian. This 
term is 

£ = a 4 \< D^ZD^ ><m Q (Z + tf) >. (21) 

Given the structure of this operator, it can be removed with a field redefinition of £ resulting 
in a multiplicative renormalization of the mass term in Eq. ([3]). This renormalization does 
not depend on the external field strength, and thus only modifies the bare pion mass. We can 
thus remove the term in renormalized perturbation theory by working with the physical 
pion mass. Thus at 0(e A ), the low-energy constants in Eq. (1201) are the only new parameters 
required. 



B. Pion in Strong Magnetic Fields 



Having written down the relevant terms of the leading and next-to-leading order chiral 
Lagrangian, we can calculate pion two-point function at 0(6 4 ). At this order, we have the 
tree-level couplings in Eq. (120]) and one-loop diagrams. The latter arise from expanding the 
leading-order Lagrangian Eq. ([3]) to quartic order in pion fields. There are three different 
diagrams generated at this order, and these are depicted in Figure [2J We will handle the 
simplest contribution first. Diagram C receives contributions from the mass term of the 
chiral Lagrangian Eq. expanded to fourth order in the pion fields. The value of the 
tadpole diagram is proportional to the self-contracted propagator, which for the charged 



3 In an electromagnetic field, where E ■ B is non-vanishing, we would additionally need to consider terms 
from the Wess-Zumino-Witten Lagrangian 33J, |34( . For a magnetic or electric field, such terms are absent 



unless other sources, such as finite density, are considered, see 



35, 
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pions we denote by T>(B). The self-contracted propagator for the neutral pion is then 
simply T>(0). In dimensional regularization, the self-contracted propagator is 



V(B)=D(x,x) = I™ dss*~ l e- sm * ^— . (22) 

v ; v ' ; (4tt) 2 J sinheSs v ; 

The result of the mass insertion vertex is a shift of the leading-order masses 

ml ± = ml — ml (l - ^[42?(B) + V(0)]j , (23) 

mlo =ml —>ml(l- ^[2V(B) + 3T>(0)}) . (24) 



The remainder of the loop contributions to the pion two-point function arise from the 
kinetic term in Eq. ([31). Expanded to fourth order, the interaction vertex consists of four 
pion fields, two with gauge covariant derivatives and two without. Diagrams A-C are related 
by gauge invariance; their contributions should be accordingly grouped in a gauge covariant 
fashion. A particularly useful simplification to note is that there are no loop contractions of 
D^ifi with 0. This is obvious for the neutral pion because the covariant derivative is simply 
a partial derivative and the loop integration must vanish by Lorentz invariance: 

J ddk ^r°- (25) 

For the charged pions, one must use the explicit form of the propagator Eq. (fT7|) to see 
that such contractions also vanish. Basically the covariant derivative of the self-contracted 
propagator ends up without coordinate dependence, hence has no four-direction in which to 
point. Thus there are only two ways the covariant derivatives can act to produce a nonzero 
result: either both act on the fields contracted to the external legs, or both act inside the 
loop. 

Consider contractions where both covariant derivatives act on the external legs. In this 
case, the remaining two pion fields are contracted to form the tadpole which by Eq. fT22|) 
is independent of the vertex location. The modification 5G(x', 0) to the two-point function 
thus has the form 

SG(x',0) = -CV(B) J d 4 xD tl D(x',x)D fl D(x,0), (26) 

where C is the overall numerical factor, and _D M acts on the x^-coordinates. We can integrate 
by parts to relocate the covariant derivative from the final state to the initial state. This 
is possible for the gauge field piece of the derivative, because in the contraction either the 
charge is zero for the neutral pion; or, for the charged pion, there must be a relative sign 
between the charges of the fields contracted to the external legs. Hence this correction to 
the two-point function has the form 

G(x',0) = D(x',0) -CV(B) J d 4 x D(x', x) (—D^D^) D(x, 0). (27) 

At this point, we could use the Green's function relation in Eq. (j!8j) . but it is best for later 
purposes to expose that we are calculating the wavefunction renormalization factor. To this 
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end, we write Eq. (I27p in an operator notation 

G(x',0) = (x'\D\0) -CV(B) J d 4 x (x'\D\x) (—D^D^) (x\D\0) 

= wc-cw-WW - ^l g i + CP ( B )(-^) O |0) ' (28) 

where, in the last line, we treat T>(B) as a perturbation. Using the identity 

1 D-^ZH + Ap 1 , (29) 



1 + AD 
we find 

G(^o)H* V 1 + cz,(VW 0) - (30) 

As an operator, D~ l = —D^D^ + m%, leading us to the final form for the correction 

G (*''»)^* \ + cvm l_ DiiDii)+ml W (3D 

The factor [1 + CT>(B)] is just the wavefunction renormalization Z. Rescaling the pion 
fields appropriately, we have 

K± — ^Z- X J 2 = vr ± (l + ^[©(S) + 2>(0)]) , (32) 



7T 



o k „o /7-1/2 _ n 



2 



^^ = ^(1 + ^(5)). (33) 



The final loop contribution arises from contractions where both covariant derivatives act 
inside the loop. In the case of the neutral pion, we have the dimensionally regulated loop 
integral 

d d k-^ = -mlV(0). (34) 

Essentially the same relation holds for the charged pions which we can demonstrate using an 
integration by parts. Using brackets between fields to denote contractions, this contribution 
from the vertex has the form 

d i x ( t>[D^D^ = J c^a:(0[-D^00]0-^[0D^]0-0[0D^]^). (35) 

Loop contractions of -D M with 0, however, vanish and the second two terms are hence zero. 
On the first term we use the Green's function relation Eq. (|T8|) and hence charged pion loop 
contributions from this term are the natural generalization of Eq. (1341) : —m\T>[B\ The 
flavor and overall normalization factors are the same as the wavefunction renormalization, 
hence that contribution is effectively doubled. 

To present the final results, we must renormalize our expressions. We choose to renormal- 
ize in terms of the physical pion mass. The net effect of this renormalization is a subtraction 
of the zero-field results. The self-contracted propagator becomes T>(B) = V(B) — V(0), 
which is free of divergences. This function can be calculated in closed form 

— , v e\B\ ( ml \ „ „ 

p ( B » = (i^ z (^)' (36) 
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where 



15 



Kx) 



ds 



sinh s 



- 1 



./■(i-io g |) +2]o A r 



l + x 



-\og2n. (37) 



Assembling the tree-level and loop contributions, the charged pion propagator at one-loop 
order has the form of D(x',x) in Eq. (j!7p . up to the replacement 



m 



cfT,^± = m ir + -p8e B (a 9 + a w ) 



(38) 



Notice that the loop contributions have exactly cancelled at this order. In net, the charged 
pion two-point function only receives corrections from neutral pion loops and these are 
unaffected by the background field. 4 The term quadratic in the external field is proportional 
to the magnetic polarizability, i.e. 



m cS (B) = m - U-kPmB 2 + 0(B 4 ). 

Comparing with this expression, we find the known next-to-leading order result [37j 

a 9 + a 10 



/ 2 m„ 



(39) 



(40) 



with the fine structure constant a/. s . = e 2 /Air. 5 Terms with B 4 or higher powers of B in the 
effective mass squared arise at two-loop order, and are hence suppressed by m^/A 2 . 

For the neutral pion, however, the situation is quite the opposite. There is no local 
contribution at next-to-leading order, and the loop diagrams do not exactly cancel. The 
neutral pion propagator at one loop maintains a Klein-Gordon form, but with an effective 
mass squared given by 



1 + 



2e\B\ 



-I 



mz 



e\B\ 



(41) 



Taking the non-relativistic limit, e\B\ m 2 , we have the perturbative expansion of neutral 
pion energy in a weak magnetic field 



E n o(p = 0) 



1 feB 
6 \ ml 



7 feB 
180 Vm! 



31 feB 



630 \mj 



(42) 



To derive this expansion, we utilized the integral representation of the function I(x), and 
series expanded the integrand. This yields only an asymptotic series because the integral 



4 This accident is additionally true to one-loop order in ST/ (3) chiral perturbation theory as well. In that 
case, the charged pion two-point function only receives loop contributions from the neutral pion, and the 
eta. 

5 Results for the electric and magnetic polarizabilities of both c harg ed and neutral pions are known at 
next-to-next-to-leading order in the chiral expansion 
recover these results in a weak field expansion. 



38. 



39. 



40. 



41|. Only at two- loop order would we 
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FIG. 3: Comparison of perturbative and non-perturbative results in constant magnetic fields. 
Plotted is the change in pion energy 5E in Eq. (|44p as a function of the expansion parameter £, 
Eq. (|45p . LO, NLO, and NNLO are used to denote the order of the perturbative approximation in 
field strength. From left to right, the three panels depict the perturbative regime, the transition 
regime, and the non-perturbative regime. 



does not converge uniformly. The first term in this asymptotic series is the energy shift of 
the neutral pion due to the magnetic polarizability, which we can read off as 



«/..< 



3(47r/) 2 m, 



(43) 



This agrees with the known one-loop result [37|, |42|, |43j . Unlike the charged pion, there are 
terms of all orders in e\B\jm\ generated in the neutral pion energy at one loop. These terms 
correspond to irreducible 7r°-multiphoton interactions. 

We can explore the behavior of the asymptotic expansion compared to the non- 
perturbative result. To do this, we consider the relative change in the neutral pion energy 
5E defined by 



5E 



(4tt/) ; 



mi 



E 7T o(p = 0) 



The exact numerical factor for the expansion parameter is arbitrary. We choose 

e\B\ 



(44) 



(45) 



so that the first two terms of the asymptotic expansion are of order unity, viz. 

^5^ 35 S 



(46) 



Because the series is asymptotic, large coefficients will automatically appear at some order 
in the expansion. We can investigate the nature of the asymptotic series by comparing the 
non-perturbative result for SE, in Eq. (144p . to the perturbative series, Eq. (1461) . This is done 
in Figure [3], where we compare the leading-order (LO), next-to-leading-order (NLO), and 
next-to-next-to-leading-order (NNLO) results for the energy shift with the non-perturbative 
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FIG. 4: Further comparison of perturbative and non-perturbative results in constant magnetic 
fields. Plotted is the change in the pion energy 5E vs. £. The dashed curve shows the non- 
perturbative result from Eq. (|4ip . The solid curves show the perturbative results to a given order, 
up to fifteenth order in the field-strength squared. The higher-order results accumulate near £ = 0. 

result. Notice we are working to a fixed order in the chiral expansion, and the expansion 
we are considering is in the strength of the external field. Three different regimes of the 
expansion are shown. In the perturbative regime, the expansion parameter £ ~ 1/4, and 
the higher-order results better approximate the non-perturbative answer. In the transition 
regime, the expansion parameter £ ~ 1/3, and adding higher-order terms is no longer ef- 
ficacious. There is a point where the LO, NLO, and NNLO results all differ by the same 
amount (in magnitude) from the non-perturbative result. This point marks the transition 
from the perturbative regime to the non-perturbative regime. The non-perturbative regime 
is typified for £ ~ 2/3. In this region, all higher-order terms push the result away from 
the non-perturbative answer. Consequently the LO result agrees best, ~ 20%, but cannot 
be improved using perturbation theory. This behavior is precisely as one expects from an 
asymptotic series: only the first few terms give a good approximation. Beyond these terms, 
additional corrections are uncontrolled. When the expansion parameter is large, we are 
necessarily restricted to just the LO term. 

Finally in Figure HJ we compare perturbative results up to fifteenth order in the expansion. 
For a given order in the expansion, there is a value for the expansion parameter £ for which 
the perturbative result diverges from the non-perturbative curve. The oscillatory nature 
of the perturbative expansion can be seen from the plot, as all odd-order approximations 
undershoot the non-perturbative result, while the even-order approximations overshoot the 
non-perturbative result. The region of £ where n th -order result gives a reasonable description 
of the exact result is bounded in magnitude roughly by 1/n. Consequently the higher-order 
curves accumulate on the plot near £ = 0, and reflect that the series expansion has strictly 
zero radius of convergence. 

C. Pion in Strong Electric Fields 

It is straightforward to repeat the calculation of the pion two-point functions at one-loop 
order in an external electric field. The diagrams which contribute are identical to those for 
the magnetic case. In general, we cannot write down the particle's energy shift in an electric 
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field from that in a magnetic field because the energy shift is not Lorentz invariant. The 
tadpole topology offers a valuable simplification, however, as the tadpole loop-factor does 
not depend on the location of the pion source and sink. For example, we need not worry 
about the relative time entering the loop integration. Because we are afforded a Lorentz 
covariant treatment to arrive at the effective mass squared, the electric and magnetic results 
are consequently the same. This only occurs to one-loop order. (Moreover this accident 
does not occur for nucleons at one loop, as we shall see below in Section [III CI ) 

Given the simple tadpole loop topology at next-to-leading order, we can thus simply write 
down the Euclidean electric field results. For the charged pion, we have 



m 7T * m cS,n± = m 7r + -p 8e £ («9 + «lo) 

= ml - j^8e 2 E 2 (a 9 + a w 



(47) 



where, in the second line, we have performed the trivial analytic continuation to Minkowski 
space. The term quadratic in the external field is proportional to the electric polarizability, 
i.e. 

m eS (E) = m-^a E E 2 + 0(E 4 ). (48) 



Comparing with this expression, we find the known next-to-leading order result [37 



a 9 + q;io 
f 2 m n 



(49) 



and the relation: a E + Pm = 0, which is a consequence of the helicity structure of the 
one-loop graphs. 

For the neutral pion, we can also write down the Euclidean space results by inspection, 
namely 



m cfr,7r° 



2eE_ fK 
(4tt/) 2 V e£ 



(50) 



but the analytic continuation to Minkowski space requires care. First let us note that the 
asymptotic expansion of i1{—ix) needed for the electric field case is just the trivial analytic 
continuation of the Euclidean space asymptotic expansion. Subtleties in the analytic contin- 
uation can only show up in non-perturbative physics. Thus we have the same perturbative 
expansion for the neutral pion energy as above under the simple replacement B — > %E (or 
£ —>«£)• The neutral pion energy consequently has the weak electric field expansion 



E 7r o(p = 0) = m n + 



eE 
ml 



7 (eE 
+ 180 \mj 



31 (eE 
+ 630 \m% 



+ 



(51) 



From this expression, we can read off the electric polarizability 

al° = ^£ 

E ~ 3(4tt/)W 

which is in agreement with the known one-loop result [37|, |42|, 



(52) 
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FIG. 5: Comparison of perturbative and non-perturbative results in constant electric fields. Plots 
are as in Figure [3] with identical regions of £. We have shown only the real part of the non- 
perturbative result. 



To investi gat e the non-perturbative result, we must perform the analytic continuation. 



As shown in [15J, the correct analytic continuation is given by 



oo 

i I(—ix) = x ^1 — 7^ — log — J — 2 tan -1 x + 

n=l 



2 tan 



-i 



x 



2n + l 



— i log(l + e 



(53) 



With this expression we can evaluate the neutral pion energy in an electric field, Eq. fl50l) . 
and compare with the weak field expansion. This is done in Figure where we compare 
LO, NLO, and NNLO results for the energy shift to the real part of the non-perturbative 
result. In each case, we plot 6E, in Eq. (jUJ), for a constant electric field. Again we are 
working to a fixed order in the chiral expansion, and the expansion we are considering is 
in the strength of the external field. We show identical ranges for the expansion parameter 
£ as for the magnetic case. Here, of course, we use £ = e\E\/ \fQrn^. In what we call the 
perturbative regime, the expansion parameter £ ~ 1/4, and the higher-order results better 
approximate the non-perturbative answer. Because all terms in the perturbative expansion 
are now positive, it must be the case that the leading-order term (and some higher-order 
terms) are less than the non-perturbative answer. In the transition regime, the expansion 
parameter £ ~ 1/3, and adding the NNLO and higher terms is no longer efficacious. This 
is where the NLO result makes its transit through the non-perturbative answer. Once 
this happens, the addition of any higher-order terms necessarily worsens the agreement. 
Increasing £ further, the LO result becomes as large as the non-perturbative answer. This 
last regime is typified for £ ~ 2/3, and is where all higher-order terms push the result away 
from the non-perturbative answer. Perturbation theory is of no avail; consequently the LO 
result agrees best. The agreement, moreover, is phenomenally good. 



6 As a warm-up for the nucleon case, we derive an equivalent integral representation for the analytic 
continuation of iX{~ix) in Appendix \C\ 
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FIG. 6: Further comparison of perturbative and non-perturbative results in constant electric fields. 
The real part of the non-perturbative energy shift is plotted along with the first fifteen perturbative 
corrections in the strength of the field squared (as in Figure 0J) . 

In Figure [6l we compare perturbative results up to fifteenth order in the expansion as we 
did above for the magnetic case. Because all terms in the perturbative expansion are positive, 
we see that for large enough £, all approximations overshoot the non-perturbative result. 
Again we are only considering the real part of this result. The higher-order approximations 
accumulate around £ = 0, as expected. Beyond the perturbative regime, the LO result 
gives a remarkably good approximation for |£| < 1. Of course this approximation cannot be 
improved. 

A feature of the non-perturbative result, is the imaginary part of the effective mass which 
we have so far neglected. This imaginary part has the form 



2 



rv / 2 \ 2e\E\m n 

^Kfluo) = - (4?r/)2 log 



1 + exp 



Ttml 



e\E\ 



(54) 



which corresponds to a width T^o = — 3m(m^ 7 ro)/ m w The full width of the neutral pion 
is a sum of partial widths 7 

oo 

TttO = ^^n, (55) 

71=1 

where the partial widths correspond to decay modes with n-pairs of charged pions 

T„ = r[vr° -?U TrV+TT-n = ("l) n+1 ^^^P (-^f ) • (56) 
Virtual fluctuations in the pion can decay to real particles in strong electric fields. 



7 This dec omposi tion follows from similar reasoning to that in the case of pair production from the vacuum, 
see, e.g., 



Ml4al4al47| 
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D. A Connection with the Chiral Condensate 



Now we demonstrate that the Gell-Mann-Oakes-Renner relation is satisfied for the pions 
in electric and magnetic fields. The magnetic result for the neutral pion was previously 



demonstrated in [12], however, that analysis was limited to the strict chiral limit, m\ = 0. 
Here we keep the pion mass non- vanishing using the power counting employed above, Eq. (J5|), 
and demonstrate the relation for both charged and neutral pions. 

In the absence of external electromagnetic fields, the Gell-Mann-Oakes-Renner rela- 
tion [48j reads 

f 2 m 2 n = 4<qq>m, (57) 

where < qq > is the quark condensate. In the chiral limit, we have f n f and < qq >^ A so 
that Eq. ([6]) is recovered. We ignore strong isospin breaking throughout so that there is only 
one condensate: < qq >=< uu >=< dd >. Now imagine turning on a constant magnetic 
field, for example. Above we have deduced the effective mass for charged and neutral pions. 
For the neutral pion, the effective mass is just the rest energy; while for the charged pion, 
this physical interpretation is not warranted because there are additionally Landau levels. 
This technicality notwithstanding, we might imagine the following generalization of Eq. (|57|) 

f 2 AB)ml cS (B)=A<qq> B m, (58) 

where we have displayed the dependence on the external field B. We can write each of the 
field- dependent quantities in terms of their zero field value and the relative difference 

<eff(£) = ml[l + 5ml(B)}, (59) 
f*(B) = f 2 Al+6f 2 AB)}> (60) 
<qq> B = <qq> [1 + 5 <qq> B }- (61) 

Taking the ratio of Eq. (1581) and Eq. (|57|) , we have 

5fl(B) + 5ml(B) = 5<qq>B. (62) 

Because each of these differences arises at next-to- leading order ~ e 2 , it is obvious that terms 
second order in differences can be dropped. To verify the conjectured relation Eq. (1581) . we 
shall merely evaluate Eq. (1621 . We have already determined the shift in the charged and 
neutral pion effective mass as a function of the magnetic field B. The relative shift of the 
condensate is known: from [15], we have 

C-^V (63) 



(4tt/) 2 \e\B 

Hence to verify the low-energy theorem Eq. (|5"8"|) . we must check that the change in pion 
decay constant compensates for the field dependence of the pion mass and chiral condensate: 

SfUB) = -^l(^), (64) 



(4tt/) 2 \e\B 
r „ 2 2e\B\^fml\ 8e 2 B 2 . 

= ~(4ir7F \e\B\) ~ T 5 "^ ' (65) 
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For the neutral pion, calculation of the decay constant in a background magnetic field is 
straightforward. The decay constant enters in the pion-to-vacuum transition matrix element 
of the axial current A a „ 

J d A x e- ip "^ (0|^(x)|tt°(0)) = VtS^ipufa. (66) 

Because the neutral pion propagator maintains a simple pole form, we can apply LSZ to 
calculate this amplitude. At this order in the chiral expansion, there are no external field- 
dependent local contributions to the axial current. This fact owes to the charge neutrality 
of the 7T°, but does not continue to hold beyond next-to-leading order. There are thus only 
two types of contributions to the external field dependence of f w o: the wavefunction renor- 
malization factor, and the loop contribution from the next-to-leading order axial current. 
The former has been determined in Eq. (133]) . The next-to- leading order axial current has 
the form 

NLO 16 

When considering the a = 3 current, recall that all the D^cf) terms must be contracted with 
the external pion, else the loop contractions vanish. Adding up the loop contributions from 
the wavefunction correction and next-to- leading order current, we arrive at 



K = — £ (D^n • t? - n a D,n ■ 7?) . (67) 



fAB) f[l-jxV{B)\ 2 



U f[l - IrP(O)] P 



[V{B) - X?(0)] + . . . , (68) 



/ 



where we have dropped higher-order terms. Squaring this result and using the definition of 
T{x) yields the desired equality in Eq. ( |64l) . 

The charged pions, by contrast, cannot be handled in this manner. One problem with 
the approach is the amputation; and, instead we focus on the coordinate-space correlation 
function C^(x), given by 

C±(x) = (0\A±(x)^(0)\0), (69) 

where it is assumed that x^ ^ to avoid contact terms. At next-to-leading order in the 
chiral expansion, there are two possible Lorentz structures into which the correlator can be 
decomposed 

C^(x)=-f 1 (B)D^D(x,0) + 2if 2 (B)F^D v D(x,0). (70) 

Terms with higher powers of the field strength tensor, such as F pv F vp D p , must come from 
local terms in the effective theory, and are hence suppressed by at least two powers of the 
chiral symmetry breaking scale. The field dependence of the coefficients fi(B) and /2(B) 
must be generated from charged pion loops. A quick power counting argument shows that 
f 2 {B) can only be generated from local interactions, hence this coefficient is independent of 
the magnetic field, /2(B) = f 2 - 

Taking the covariant four- divergence of the axial correlator, we have 

D»C±(x) = -[f x {B)ml + f 2 QB 2 }D(x, 0). (71) 

To arrive at this expression, we have used the fact that the propagator is a Green's function 
Eq. (TT81) . and the relation 2F pu D fl D l/ = iQF^ v F pv . In this form, we can amputate the 
correlation function in coordinate space 

(01^^(0)) = Jd*x(-D^ + rn^)D^(x). (72) 
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From this amplitude arises a natural generalization of the pion decay constant 



(0\D,A±\n±(0)) = -mlU(B). (73) 

In the chiral limit, the axial current is conserved d^A a ^{x) = 0. Thus the decay constant 
defined in Eq. (173]) will be singular in the chiral limit because the covariant divergence does 
not vanish. We are forced to use the covariant divergence, however, in order to utilize LSZ. 
Notice that the amplitude itself is finite in the chiral limit. The singularity in f n ±(B), 
moreover, must disappear when the external field is turned off. 

Having generalized the defining relation for the pion decay constant, we can now calculate 
f n ±(B) to one-loop order. The wavefunction correction has been obtained above, Eq. (1321) . 
Additionally one-loop diagrams arise from the next-to- leading order axial current in Eq. (1671) . 
Combining the results from loop diagrams, we can identify the factor fi(B) in Eq. (1701) . 
namely 

/i(B) = /|l--i[P(B)+2?(0)]}. (74) 

The coefficient f'2 arises from a local contribution to the axial current contained in the 
next-to-leading order Lagrangian. The relevant external field dependent terms are just those 
in Eq. (120]) . but generalized to include both left- and right-handed sources 

£ = -ia 9 < L^D^D V Y) + R^D^DJZ > -a w < L^ZR^tf > . (75) 

The left-handed field strength, L^ u , is given by L^ u = d^L u — d v L^ + i[L^, L u ], and similarly 
for the right-handed field strength, R^ u . The left- and right-handed four- vector potentials, 
L M and R^, are given by: L M = QeA^ + T a A a , and R^ = QeA^ — r a A a . The action of the 
chirally covariant derivative is specified by, D^Y, = + — iHR\ Evaluating the single 
pion terms of the axial current derived from Eq. (175]) . we find the tree-level contribution 

SA* = -j8i(a 9 + a 1Q )QeF IJa ,D v 'K ± , (76) 

with SA^ = 0. From this expression, we can calculate the axial correlator and identify f'2 
in Eq. (1701) . We find the value f 2 = — 4(a 9 + Oi W )Qe/ f. Combining the tree-level and loop 
contributions, we deduce the charged pion decay constant using Eq. (1731) 

f ±(B) 1 4e 2 B 2 

= 1 " ±[V{B) - V(0)) - ^4(a 9 + a w ). (77) 

Upon squaring, we find that Eq. (165]) is indeed satisfied. 

We have shown that the low-energy theorem for the neutral pion in a magnetic field 12| 



extends beyond the chiral limit into the regime where m 2 ~ e\B\. Furthermore we have 
generalized the definition of the charged pion decay constant for a background field. This 
generalization enabled us to demonstrate that the Gell-Mann-Oakes-Renner relation holds 
additionally for charged pions in magnetic fields, Eq. (1581) . There is no subtlety in repeating 
the calculation for the case of a Euclidean electric field. The effective mass was deduced 
above in Section Hi CI The chiral condensate is Lorentz invariant, and the determination of 
the decay constant was Lorentz covariant. Thus we arrive at the £-field results for < qq >s 
and f n (£) merely by substituting B —>■ S. The final step to deriving the low-energy theorem 
in an electric field is the analytic continuation. The tree-level contributions are trivial to 
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continue, while the one-loop expression requires Eq. (1531) . After analytic continuation, we 
arrive at 

fl{E)mlJE)=A<qq> E m, (78) 



for 7r = 7T°, and tt ± . 



III. NUCLEON IN STRONG FIELDS 
A. Nucleon Chiral Lagrangian 

To describe the physics of nucleons in strong electric and magnetic fields, we use the 
nucleon chiral Lagrangian. A well known complication arises in trying to describe the chiral 
dynamics of the nucleon. In the chiral limit, the nucleon has a non- vanishing mass, M^, 
which is on the order of the chiral symmetry breaking scale, ~ A x . Consequently there 
is no power counting to order the tower of operators needed to renormalize the theory of 
nucleons and pions, e.g. terms with any number of derivatives acting on the nucleon field 
are O(e ). The solution to this complication is to treat the nucleon non-relativistically. A 
particularly clean implementation of the non-relativistic expansion for nucleons is known as 



heavy baryon chiral perturbation theory [49j, |50| , which borrowed earlier developments from 
heavy quark effective theory [5l|, 52| . 

Consider a fermion field B of mass Mg, with Mg = M^ + A. We choose this consideration 
because we will include both nucleons (for which we can choose A = 0), and deltas (for which 
consequently A ^ 0) into chiral perturbation theory. The free Lagrangian for B in Euclidean 
space appears as 

C = B(p + M B )B. (79) 

The mass terms spoils the power counting when pion interactions are included. To remove 
this term, we expand about the non-relativistic solution to the classical equations of motion 

B v (x) = V + e- iMNV ^B{x). (80) 

Here is the covariant four- velocity. In the rest frame, = (0, 0, 0, i), so that the velocity 
satisfies the Euclidean space relation v^v^ = —1. The conjugate field B v satisfies the conju- 
gate relation without ty, complex conjugated. In Euclidean space, we are free to treat the 
field and its conjugate as independent variables. The projector V+ is given in an arbitrary 
frame by V+ = |(1 — if). Inserting the field redefinition, Eq. ( IHOI) . into the free Lagrangian, 
Eq. (J7HD, we find __ 

C = B v {-iv^d fl + A)B v . (81) 

Having phased away the large mass Mn, we are free to treat the residual mass A as a 
small parameter. Additionally with the Fourier modes, of the original B field written as 

= Mv^ + k^, we see the residual momentum k^ is produced by derivatives acting on B v . 
For slowly moving baryons, all components of are necessarily small. 

Before writing down the heavy nucleon chiral Lagrangian, let us derive the free propagator 
for the B field in the static limit. The heavy baryon two-point function, Dg(x, 0), has the 
form 

D B (x,0) = (0|i3(x)S(0)|0) = e lM «(0|P + H l ,(x)^(0)P + |0) 

TTV.-V 1 a = V + 5(x)9(x,)e-^ +A ^, (82) 
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where, after the first line, we have used the rest frame. This propagator is the essential new 
ingredient needed to perform chiral loop calculations for the nucleon. 

To describe the nucleon, we use a two component isospinor iVj, given by N = (p,n) T . 
Additionally we include the nearest lying resonances, the deltas. These give rise to generally 
important virtual corrections to nucleon observables in chiral dynamics, because the mass 
splitting A between the nucleons and deltas is on the order of the pion mass. Furthermore 
the pion-nucleon-delta axial coupling is not small. The deltas we describe using a completely 
symmetric flavor tensor T^, where: T m = A ++ , Tn 2 = A + />/3, T 122 = A /a/3, T 222 = A - , 
etc. To keep a consistent power counting, we must expand the N and T fields about their 
classical equations of motion. We shall leave the velocity subscript as implicit. Because the 
nucleons and deltas interact, we cannot phase away both chiral limit masses and Mt- 
Instead we phase away Mjv from both multiplets. This leaves a residual mass A = Mt — Mjy 
for the deltas. Thus we have the following additions to the power counting in Eq. <^ 

A ~e. (83) 



M N M N 

We will not need to distinguish between the heavy baryon mass, and the chiral symmetry 
breaking scale, i.e. we treat M N ~ A x . 

The relevant terms of the heavy nucleon and delta chiral Lagrangian at 0(e) are 

C = A7(-iV^+^HV^+ A )^ ( 84 ) 

The vector is the covariant spin operator satisfying S^S^ = 3/4. Electromagnetism is 
coupled to the Largrangian through the axial A^ and vector V M fields of pions given by 

V, = ieA^Q + -L (0^0 -D^)..., (85) 
1 

T 

where higher-order terms have been dropped. The chirally covariant derivative X> M has the 
action 

(2VV) . = d^N t + (V M ), l ' N v + tr (V M ) N u (87) 

on the nucleon field, and 



Afj, = -D^ +-.., (86) 



{V„T) m = d^Tijk + (V M ), l 'T V3k + (V M ) . fTwu + (V M ), k 'T ijk ,, (88) 



'j 

on the delta. 



B. Nucleon in Strong Magnetic Fields 

Having outlined the chiral Lagrangian for the nucleon, we now determine the energy of 
a nucleon in an external magnetic field to 0(e 2 ) by computing corrections to the nucleon 
two-point function. Corrections to the nucleon's energy in an external field arise from pion 
loops. As in the meson sector, these contributions are non-perturbative in the field strength, 
e\B\/m% ~ e°. There are additionally local interactions, but the coefficients of these terms 
are phenomenologically known. Thus we will be able to make quantitative predictions for 
nucleons in strong fields. 
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FIG. 7: Diagrams for the nucleoli two-point function in an external field. The straight lines denote 
baryons, while the dashed lines represent the background field propagator of the pion. Wiggly lines 
terminating in crosses represent couplings of the interaction vertices to the background field. The 
internal baryons lines can be both nucleons and deltas. 

In calculating the nucleon two-point function to 0(e 2 ), there are both local contributions 
and loop contributions. In the effective field theory, the local interactions arise from magnetic 
moment operators. For the nucleons, there are isovector and isoscalar operators having the 
form [H, 54 

ip 

C = ~^r F ^ {t*o N [ S » S U }N + ^N[S„ S u }t 3 N) . (89) 

There is only one additional operator needed to this order: the nucleon-delta magnetic dipole 
transition operator [53( 

£ = sVSfir-^-F^ [(NS,QT U ) + (T^QSvN)} . (90) 

AIVIn 

In the absence of other corrections, these terms lead to a shift in the nucleon energy of the 
form 

The magnetic moments denoted by /xjv are isospin dependent, /i p = f-io+f-ii, and fi n = fiQ— 
The nucleon magnetic moment term scales as a relative 0(e) correction to the mass Mat, 
and the magnetic transition term scales as 0(e 2 ). The isovector magnetic moment n\ also 
serves as a counterterm to renormalize the 0(e 2 ) loop contributions. Technically one needs 
additional operators having exactly the same structure as those in Eq. (I8"9l . but multiplied 
by the chiral singlet parameter A/M^r. In addition to renormalizing contributions to the 
nucleon mass, we will additionally renormalize the magnetic moment contributions so that 
Hi takes its physical value. 
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There is one further local operator that leads to a nucleoli energy shift in an exter- 
nal magnetic field. This correction has a coefficient exactly fixed by gauge invariance and 
reparametrization invariance [55j|. The latter arises as a convenient way to impose Lorentz 
invariance order-by-order in the heavy baryon expansion. Applying reparametrization in- 
variance on the leading kinetic term for the nucleon, we generate the kinetic energy operator 



where (T>±)n = T)^ + v^vJDv). This operator leads to corrections to the nucleon energy 
at 0(e 2 ), thus we do not need the chiral structure in the derivative X> M . Furthermore, for 
the neutron, we completely eliminate this contribution by keeping the neutron at rest. For 
the proton, however, the gauged kinetic energy leads to Landau levels. For zero transverse 
momentum, the proton Landau levels are given by 



where n is an integer. Although these energies scale as O(e), large values of n will upset 
the power counting, so we are necessarily restricted to only the lowest levels. Relativistic 
corrections to this result give rise to an energy shift at 0(e 3 ) in the power counting. 

The remaining contributions to the nucleon energy are non-analytic in nature, and stem 
from pion loop diagrams generated from chiral dynamics. The interactions between baryons 
and pions are contained in Eq. ( 1841 . Using the vertices appearing in this Lagrangian, 
we generate the diagrams shown in Figure [7] for the nucleon two-point function. Most 
importantly there are no charge couplings of the photon to the nucleon or delta. These 
have been eliminated from Eq. (1841) by the gauge condition v^A^ = 0. Consequently the 
background field only couples to the pions, or to the interaction vertices. 

Of the diagrams shown in the figure, five are trivial to evaluate. Diagrams A\ and A2 
arise from expanding out the vector field of pions to include two pion terms. The partial 
derivative piece generates diagram A±, while the gauge piece of the derivative generates 
A 2 . Both of these diagrams vanish because the vector field is dotted into the four- velocity 
v^. In the case of A\, the derivative acting inside the loop picks off a good component of 
momentum, making the momentum integrand odd; in the case of A 2 , we have v^A^ = 0. 
While the remaining diagrams B1-B3, and C are all connected by gauge invariance, we 
find it simpler to evaluate each separately. A useful simplification to note stems from the 
Dirac-delta function in the heavy nucleon propagator, Eq. f|82|) . We can choose the nucleon 
source to be located at position x = 0. Because the nucleon remains static, the Dirac-delta 
function enforces the intermediate and final states also to be at x = 0. For diagrams B1-B3 
coupling of the external field to the vertices introduces an explicit power of the potential A^, 
which depends on the spatial coordinates. These diagrams thus vanish upon imposing the 
Dirac-delta function on position. The only non-vanishing loop contribution to the nucleon 
two-point function in an external magnetic field is from diagram C. As the evaluation is 
rather technical, we relegate the calculation and subsequent renormalization to Appendix [Al 

The complete expression for the nucleon energy in a background magnetic field up to 
0(e 2 ) is given by 




(93) 



E N (B) =M N + 
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2M N 




+ 5E 2 (B 2 ). 



(94) 
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Here r 3 is the usual isospin matrix, r 3 = diag(l, —1), and we have specialized to the case of 
the lowest Landau level. The non-analytic functions 5Ei(B 2 ) and 5E2(B 2 ) arise from pion 
loop graphs, and are given by 



5Ei(B 2 



5E 2 {B 2 ) 
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(96) 



Notice we have renormalized the function SEi(B 2 ) so that 5Ei(0) = 0. In this renormal- 
ization scheme, /i^ is the physical nucleon magnetic moment. Thus we have absorbed the 
well-known one-loop chiral corrections to magnetic moments that are linear in m n into the 
coefficient //jv- This is possible because we are considering the external field dependence and 
not the pion mass dependence of the nucleon energy. 

In weak fields, we can perturbatively expand these results to arrive at the asymptotic 
series 



5E 1 (B 2 ) = Y^dEf^B^, and 5E 2 {B 2 ) = ]T 5Ef n) B 2n . 



(97) 



n=l 



n=l 



The first term in the expansion of 5Ei(B 2 ), which gives a contribution to the nucleon energy 
at 0(B 3 ), is 
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In the chiral limit, this coefficient is singular and consequently highly constrained by the 
parameters in the chiral Lagrangian. In general, the chiral limit behavior of terms in the 
asymptotic expansion are 
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The first term in the expansion of 5E 2 (B 2 ) gives a contribution to the nucleon energy at 
0{B 2 ) of the form 
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This energy shift is the chiral loop contribution to the nucleon magnetic polarizability, 
namely /?j^° p = —^dE^K The remaining contribution to (3m comes from the transition 
moment term proportional to fi^ i n Eq. (1941) . The derived value for (5m agrees with the known 
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FIG. 8: Comparison of perturbative and non-perturbative results for relative shift in unpolarized 
neutron energy in constant magnetic fields. Plotted is the non-perturbative result for < AE n > 
in Eq. (|1U2|) as well as the LO and NLO weak-field approximations to this result. The expansion 
parameter £ = e\B\/m%. 



chiral perturbation theory expression 56|, |57|, |58|. The next correction in the expansion of 
5E2(B 2 ), gives a shift at 0(B A ) to the nucleon energy of the form 
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From these first few terms in the weak-field series expansion, we can compare the asymp- 
totic expansion with the non-perturbative energy shift. The energy of an unpolarized neu- 
tron has an expansion in only even powers of the magnetic field. We consider the relative 
shift in the unpolarized neutron energy 



< AE n >- 
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(102) 



The LO expansion of < AE n > is due to the polarizability [from the transition moment and 
Eq. fllOOp ]. while at NLO we add to the polarizability the fourth-order term in Eq. (11011) . 
We can assess the behavior of the asymptotic series using phenomenological input for the 
low-energy constants. For the nucleon-delta mass splitting, we take A = 0.29 GeV. The 
nucleon axial couplings is well known from experiment, qa = 1.25. The nucleon-delta axial 

The dipole transition 




In Figure El 



coupling, (Jan = 1.5, we estimate from the width r(A — > Nn) [59 1 
coefficient is also known from A — > A^7, giving the value ht = 1.4 
we plot the non-perturbative shift < AE n > along with the LO and NLO perturbative 
results. The expansion parameter £ is taken to be £ = e|-B|/m 2 . The two graphs contrast 
the perturbative regime, ^ C 1, where the higher-order correction gives better agreement 
with the non-perturbative result, and the non-perturbative regime, £ > 1, where LO is the 
best available approximation but cannot be improved. 
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FIG. 9: Comparison of perturbative and non-perturbative results for the relative shift in the 
spin averaged, isospin averaged nucleon energy in constant magnetic fields. Plotted is the non- 
perturbative result for < AEj > in Eq. (|1U3|) as well as the LO and NLO weak-field approximations 
to this result. The expansion parameter £ = e\B\/m 2 . 



The polarized isovector matrix element of the nucleon has an expansion in only odd 
powers of the magnetic field. Taking the isospin average, and average of spin-up minus 
spin-down contributions, we find the relative shift in the energy 

The LO term in the expansion is just the linear dependence on the magnetic field due to the 
isovector magnetic moment, /xi = 2.35. The NLO term adds on the B 2 correction given in 
Eq. (198]) . In Figure [9j we plot the non-perturbative shift < AEj > along with the LO and 
NLO perturbative expansions. The expansion parameter £ is again taken to be £ = e\B\/m 2 , 
and we depict the behavior of < AEj > in the perturbative and non-perturbative regimes. 

The proton-neutron splitting is an interesting quantity to consider in strong magnetic 
fields. Taking the difference in ground state energies of the proton and neutron, we have 

AE pn = E P (B) - E n (B) 

= M p - M n - f fio - I ] . (104) 



M N V 2 



Notice that both loop contributions, SEi(B 2 ) and 5E 2 {B 2 ), cancel out leaving the estimate 



of [63j. The first chiral loop contributions to AE occur from the spin- dependent isoscalar 
terms, and spin-independent isovector terms both present in the 0(e 3 ) expression for the 
energy, whereas here we work to 0(e 2 ). Using the value for the isoscalar magnetic moment, 
/io = 0.44, we find that the nucleon is stable to weak decay, AE pn = 0, for £ = 1.0. 
This corresponds to a magnetic field strength of B = 3.3 x 10 14 T. Proton beta decay to 
neutron, positron, and neutrino occurs when AE pn m m e , or when £ = 1.5 corresponding to 
B = 4.9 x 10 T. It is interesting to note that these values of the field strength are exactly 
in the regime where the power counting in Eqs. (jSJ) and fl83|) applies. Corrections to these 
estimates from 0(e 3 ) terms are proportional to m^/M^ ~ 15%. 
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C. Nucleon in Strong Electric Fields 



We now determine the one-loop corrections to the nucleon two-point function in an ex- 
ternal electric field. We begin first by considering the Euclidean space electric field specified 
by the gauge potential in Eq. (jSJ). After deriving the two-point function, we address the 
analytic continuation to Minkowski space. 

The first point to note is that the nucleon energy in an external Euclidean electric field 
is not directly related to that in a magnetic field. This should be obvious from Eq. (1941) . 
for example, there is no nucleon electric dipole moment. At a deeper level, while the gauge 
fields only differ by a naming of the coordinates, the nucleon correlation functions require 
us to identify single particle contributions. In the presence of background magnetic and 
electric fields, the pion fluctuations cohere differently in time to form a nucleon, and this is 
borne in upon LSZ reduction. We must thus perform a separate calculation to arrive at the 
electric field result. 

There are no local contributions to the nucleon two-point function in an electric field until 
0(e 3 ). The loop diagrams contributing to the electric-field dependence of the nucleon two- 
point function to 0(e 2 ) have been depicted above in Figured Diagrams A\ and A2 vanish; 
A2 for precisely the same reason as with the magnetic case, v^A^ = in our gauge. The 
derivative, v^d^, acting inside the loop in Ai does not bring down a power of good momentum 
but the Gaufiian integral is odd after a shift of integration variable. Diagrams B2 and B% 
are also simple to evaluate. Insertion of the gauge potential produces the time at the vertex. 
Upon amputation of the external propagators, the vertex times are set equal to the times of 
the source and sink (and only the relative time enters the meson propagator). Locating the 
source at t — forces these diagrams to vanish. The remaining two diagrams, Bi and C, 
yield non-vanishing contributions to the nucleon two-point function. Their evaluation and 
renormalization have been relegated to Appendix IB1 

The complete expression for the nucleon energy 8 in a Euclidean electric field £ to 0(e 2 ) 
is given by 

E N {£) = M N + g 2 A 5E(0,£ 2 ) + ^g 2 AN 5E(A,£ 2 ). (105) 

The function 5E(A,£ 2 ) represents the non-analytic contributions from loop diagrams, and 
has the explicit form 
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In evaluating the heavy nucleon two-point function in an electric field at this order, we can still speak 
of the nucleon energy. For the proton, in strong fields we must ultimately sum the Born couplings to the 
charge. In that case, we can still speak of an effective energy, as was the case with the charged pion above. 
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Series expanding the integrand in powers of the field 8, we arrive at the asymptotic series 



5E{A,£ 2 ) = ^5E {2n) {A)£ 2n . 



(107) 
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From this expansion, we can read off the terms in the asymptotic expansion in Minkowski 
space 

oo 

n2n 



6E(A, —E 2 ) = ^(-l) n 5£ (2n) (A) E 2 



(10£ 



n=l 



Comparing with the general form of the energy in a weak electric field, 
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(109) 



we can identify the electric polarizability a E as the first term in the asymptotic expansion, 
namely 

1 \g 2 A SE^(0) + 8 -g 2 AN 6E^(A) 



9' 



110) 



Deriving the first coefficient in the asymptotic series from Eq. (11061) . we find 
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from which it follows that SE^ 2 \0) = 5/(967r/ 2 m 7r ). From this term in the ex pan sion, we 
deduce the known one-loop value for the nucleon electric polarizability (Hfjl, 57,113]. Notice 
this first coefficient in the weak field expansion is singular in the chiral limit. In general, 
the behavior of the coefficients near the chiral limit is: <5£ (2n) (A) ~ m\ . Working to one 
order higher in the series, we find the coefficient of the fourth-order term 
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from which follows the value ^^(O) = — 5/(2567r/ 2 m^.). 

From these first few terms in the perturbative expansion in field strength, we can contrast 
with the non-perturbative result in Eq. (I105p . This comparison with the LO result from 
Eq. (flTTi) . and NLO result from both Eqs. (ITTTjl and (111211 is shown in Figure HD1 We 
show two distinct regions for the expansion parameter £ = e\E\/m 2 K . When £ < 1/3, the 
perturbative expansion is under control. For £ > 1, the expansion has broken down, and 
the LO result give the best approximation to the non-perturbative answer. 

It remains to consider the analytic continuation necessary to arrive at the nucleon energy 
in an electric field E. The perturbative expansion has been trivially continued in Eq. (11081) . 
Unlike the pion case, however, the non-perturbative result for E^{8) has a rather compli- 
cated structure in the complex plane. The analytic continuation can be performed using a 
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FIG. 10: Comparison of perturbative and non-perturbative results for the nucleon energy in con- 



stant Euclidean electric fields. Plotted is the non-perturbative result for AEn 



E N (£)-M N 
M N 



, with 



En(£) given in Eq. (|105p . as well as the LO and NLO weak-field approximations to this result. 
The expansion parameter £ = e\£\/m 2 . 

trick, as detailed in Appendix O As the resulting expressions are rather cumbersome, we 
restrict our attention to the case A = 0. In that case, we have from Eq. ( 1C9I) 
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coshs\/coths 

From the imaginary part of 5E(A, —E 2 ), one can deduce the total nucleon decay- width in 
an electric field. Here p (n) decays to n + n (n^p) plus any number of 7r + 7r~ pairs. 

In Figure HH we compare results in Minkowski space for the shift in nucleon energy as 
a function of the expansion parameter £ = e\E\/m 2 . We plot the real part of the non- 
perturbative result for the relative shift in energy AE = - M '— — . For simplicity, we have 
artificially set the nucleon-delta axial coupling to zero, = 0, to exclude contributions 
from virtual deltas. With this simplification, we can easily derive higher-order terms in the 
perturbative expansion, Eq. (11081) . and have done so up to fifteenth order, n = 15, in the 
figure. The higher-order results accumulate near £ = 0, which is by now familiar behavior. 
Additionally the LO shift in nucleon energy proportional to E 2 does well to describe the 
non-perturbative result not only qualitatively, but sometimes quantitatively, beyond the 
perturbative regime. 
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FIG. 11: Further comparison of perturbative and non-perturbative results for the nucleon in con- 
stant electric fields. The real part of the non-perturbative energy shift (for g&x = 0) is plotted 
along with the first fifteen perturbative corrections in the strength of the field squared. The ex- 
pansion parameter £ = e\E\/m%. 

IV. CONCLUSION 

Above we have investigated the effects of strong electric and magnetic fields on hadrons 
using chiral perturbation theory. We worked with a power counting that treats the strength 
of the external field non-perturbatively, while maintaining the perturbative chiral expansion. 
The field strengths meeting this criterion are B < 10 15 T for magnetic fields, and E < 
10 24 V/m for electric fields. 

For fields in this regime, we derived the effective mass of the charged and neutral pions. 
Fortuitous cancelations at one loop lead to only a quadratic shift of the charged pion effective 
mass. This shift encompasses the electromagnetic polarizability. For the neutral pion, loop 
graphs give rise to an effective mass with non-analytic dependence on the external field. The 
neutral pion polarizability is the first field-dependent term in the asymptotic expansion of 
the non-perturbative result. The calculation of the pion effective mass was easily carried 
out for a Euclidean electric field, owing to simplifying topological features of the one-loop 
graphs. These results were then analytically continued to Minkowski space. Asymptotic 
expansions were explored in the weak field limit, showing remarkable agreement between 
the non-perturbative result and the leading-order contribution well beyond the range where 
the perturbative expansion in field strength breaks down. This agreement, however, cannot 
be improved upon by higher-order terms in the series. 

We considered the effects of strong fields on the nucleon by using heavy baryon chiral 
perturbation theory, and including the nearby delta resonances. In magnetic fields, we 
derived the shift of the nucleon energy accounting for both spin and isospin dependence. 
Expressions for polarized and unpolarized energy shifts that are non-perturbative in the 
field strength were derived, and compared to their perturbative weak-field limits. Chiral 
analysis confirmed that for sufficiently strong magnetic fields the proton becomes unstable 
to beta-decay. The energy difference between ground-state nucleons was demonstrated to 
be largely insensitive to non-perturbative effects from the external field. The first chiral 
corrections to the ground-state energy difference occur at one order higher than we worked. 
It would be interesting to determine such corrections. Finally we considered the effects of a 
strong electric field on the nucleon. An intricate calculation using the coordinate-space LSZ 
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reduction allowed us to calculate the field-dependent shift in nucleon energy. The result 
exhibits considerably complicated analytic structure. We compared the non-perturbative 
Euclidean space result to the weak field expansion. Using a novel trick, we were able to 
analytically continue the nucleon result from a Euclidean electric field to Minkowski space. 
The real and imaginary contributions to the nucleon energy were expressed in terms of well- 
defined one-dimensional integrals, albeit still fairly complicated. Simplifying to the case 
without delta loop contributions, we then compared the real part of the energy shift to the 
perturbative expansion in Minkowski space. Similar to the pion case, the leading-order term 
in the perturbative expansion in field strength agrees well with the non-peturbative result 
when perturbation theory does not apply. Again improvement is not possible: higher-order 
terms in the strength of the field necessarily spoil the agreement. 

A practical application of our work concerns background field measurements of magnetic 
moments and electromagnetic polarizabilities using lattice QCD simulations. In the back- 
ground field approach, one simulates QCD with classical electromagnetic fields. The hadron 
spectrum is measured at various strengths of the external field. The measured energy shifts 
are then used to deduce electromagnetic properties. On a torus, constant fields must satisfy 
quantization conditions which accordingly restrict the available values of field strengths. On 
current-size lattices, the allowed field strengths are rather large. Furthermore as pion masses 
are brought down to the physical point, the expansion parameter, e\B\/m% or e\£\/m%, will 
approach or even exceed unity. It is thus important to understand non-perturbative effects 
the external field has on QCD bound states. Knowledge of higher-order terms in the field 
strength, and their behavior in the chiral limit, will additionally aid in a cleaner extraction 
of magnetic moments and electromagnetic polarizabilities from lattice data. 

Finally based on our study, one can make some rather speculative observations about 
perturbative expansions. There are well-known cases where perturbation theory is ap- 
plied in regimes of questionable validity. For example, 577(3) chiral perturbation theory 
for baryons 9 has an expansion parameter of questionable smallness, e ~ m^jMs ~ 1/2. 
The leading SU (3) corrections to the octet baryon masses lead to the Gell-Mann-Okubo 
formula, which experimentally holds to less than one percent. Could this auspicious result 
be the consequence, not of a controlled expansion about the SU(3) chiral limit, but of a 
perturbative expansion that has already broken down? In considering strong external elec- 
tric and magnetic fields, we have seen a few examples of leading-order terms in remarkable 
agreement with non-perturbative results. This agreement, of course, must be observable de- 
pendent because such behavior is not characteristic of all asymptotic expansions. If SU(3) 
chiral perturbation theory has broken down at the physical strange quark mass, we might 
expect some leading-order results to work well, while others would show uncontrolled be- 
havior. Phenomenologically it is rather difficult to assess the convergence of the three-flavor 
expansion, but with lattice QCD data, we may soon have new insight. 

Lattice QCD extrapolations themselves provide another example where the applicability 
of a perturbative expansion is questionable: how light must the pion mass be in order to 
extrapolate using chiral perturbation theory to the physical point? Interesting lattice results 



have been obtained for ttit, tiK, and KK scattering lengths [64], [65], [66|, |67|. The lattice- 
determined scattering lengths at values of the pion mass ranging between two and three times 
the physical value appear to be in remarkable agreement with current-algebra predictions 



made long ago by Weinberg [68fl. Viewed from a modern standpoint, such predictions are the 



9 The meson sector is not necessarily perturbative either: the expansion parameter e 2 ~ m^/A 2 , 1/4. 
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leading-order contributions from tree-level graphs in chiral perturbation theory. We might 
speculate that this tenacious agreement with leading-order chiral perturbation theory arises, 
not from surprisingly small next-to-leading order corrections, but from the pion masses being 
too large for the perturbative expansion to be reliable. Ultimately lattice QCD simulations 
nearing the physical point will shed light on the regime of validity of chiral perturbation 
theory. Nonetheless, the study of hadrons in strong electric and magnetic fields gives us an 
exactly soluble framework for studying non-perturbative effects, and gaining insight in to 
the behavior of perturbation theory. 
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APPENDIX A: TECHNICAL DETAILS: NUCLEON IN MAGNETIC FIELDS 



In the main text, we argued that the only contributing diagram to the nucleon propagator 
in an external magnetic field is that of topology C in Figure [71 Here we detail the calculation 
of this diagram, and its renormalization. 

Let G(x, 0) denote the full nucleon two-point function, and Dg(x, 0) the free heavy baryon 
propagator. The contribution SG(x, 0) from the sunset diagram with an internal nucleon 
has the form 



(Al) 



5G(x,0) =C / <ty / d*zD N (x,y)S- d y D(y,z)S- d z D N (y,z)D N (z,0), 



where C is a flavor dependent coefficient. Recall D(y,z) is the meson propagator in the 
background field, Eq. ffTTl) . We have used subscripts on partial derivatives to denote which 
variables they are derivatives with respect to. In the rest frame, the spin operators are 
purely spatial, and we thus need only calculate spatial gradients of the meson propagator. 
The spin algebra reduces to spin diagonal and non-diagonal terms via the identity 
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A similar decomposition additionally holds for the diagram with an intermediate-state delta. 
The spin factors, however, are different owing to the spin-3/2 propagator Py, which satisfies 



2 2 . 



(A3) 



The heavy baryon propagators simplify the loop integral because the spatial coordinates, 
y and z, are forced to zero (the location of the source). Focus first on the spin-diagonal 
term. The action of the partial derivatives on the meson propagator produces 
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where a = tanh ^QBs, and Q is the charge of the loop-meson. Eq. (1A4j) is an even 
function of Q, so we only need to distinguish between charged and neutral pion contributions. 
Evaluation of the momentum integrals, with Q set to unity for simplicity, yields 



V y -V z D(y,z) 



y =z=o (An) 2 J sinh£>s 



sB coth Bs + - 



(A5) 



in d = A — 2e dimensions. This result depends on the relative time T between the vertices, 
given by T = y 4 — z±. The intermediate-state baryon propagator also depends on T. Upon 
amputation of the diagram (which is possible in coordinate space due to the simple form 
of the heavy baryon propagator), this relative time turns into the time difference between 
source and sink. To find the self-energy correction, we project the nucleon energy to zero 
by integrating over the relative time. The required relative-time integral has the form 



dTe-^-^ = Erfc (A6) 

in the case of a delta intermediate state. Erfc (a;) is the complement of the standard error 
function. The nucleon result follows from taking A — > 0, for which Erfc(0) = 1. To complete 
the evaluation of diagram C, we assemble the result of the loop integration with the spin and 
flavor factors. The zero-field result gives the perturbative correction to the nucleon mass. 
Renormalizing to the physical nucleon mass amounts to a simple zero field subtraction. The 
result was given previously as 5E 2 (B 2 ) above in Eq. fl96|) . 

The non-diagonal spin term is similarly evaluated. The meson propagator evaluated at 
zero spatial position is given by 

V y . VMV, ») - r ?l^ e -Un-^X (A 7) 



y=z=o (An) 2 J sinh QBs 

This contribution is odd with respect to the meson charge and we must keep track of 
each pion contribution separately. Amputation and projection onto zero energy proceeds 
just as with the spin-diagonal case; the relative time integration is identical. Lastly this 
contribution needs to be renormalized. The linear B term produces the one-loop correction 
to the magnetic moment, which is subtracted by renormalizing to the physical magnetic 
moment. The resulting shift of the nucleon energy, we have denoted by 5E\(B 2 ) and is 
given in Eq. fl95|) above. 



APPENDIX B: TECHNICAL DETAILS: NUCLEON IN ELECTRIC FIELDS 

In the main text, we found that the only contributing diagrams to the nucleon propagator 
in an external electric field are Bi and C shown in Figure [7J Here we detail the calculation 
of these diagrams, and their renormalization. 

First we handle diagram B\. The spin factors are easy to handle, there is only a spin 
diagonal term. The combination of vertices produces S 3 S- V. The spatial gradient produces 
factors of good momentum when acting on the meson propagator in Eq. ffl9l) . Thus when 
integrating over the momenta, only V3 can be non-vanishing because the integrand is not a 
simple Gaufiian in fc 3 momentum. Hence we are lead to the spin diagonal parts of Eqs. (IA2j) 
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and (1A3I) . By virtue of the static approximation, the propagator is evaluated at initial and 
final coordinate locations situated at the origin. The action of the derivative on the meson 
propagator produces 
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We then multiply by the heavy baryon propagators, amputate the external legs, and integrate 
over the relative time between source and sink. The latter requires the integral 
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This contribution is finite and vanishes in the zero field limit, hence there is no renormal- 
ization required. 

To evaluate diagram C, we need to consider the action of two factors of S ■ V on the 
meson propagator. Because the k\ and ki integrals are simple Gaufiians, Vi^Vj acting in 
the loop is proportional to a Kronecker-delta, for j = 1, 2, and 3. Thus there are again only 
spin-diagonal contributions from Eqs. (IA2I) and flA3j) . The relevant contribution from the 
meson loop is 
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with T>(£) defined in Eq. fIBlj) . To this expression, we append the propagators for the 
heavy baryons, amputate, and finally integrate over the time difference. The relative time 
integration requires Eqs. (IA6I) and (IB2I) . The resulting expression diverges only in the 
zero field limit. Subtracting off the divergence amounts to renormalizing the nucleon mass. 
Performing the subtraction and combining the results of diagrams B\ and C with the relevant 
flavor factors, we arrive at Eq. ( 11051) given in the main text. 



APPENDIX C: ANALYTIC CONTINUATION 

In calculating hadronic two-point functions in external electric fields, we utilized a Eu- 
clidean electric field £. To arrive at the Minkowski space result requires an analytic contin- 
uation in the field strength, £ — > iE. For the pion case, the requisite continuation rather 
fortuitously exists in closed form and is given by Eq. (jMD, see flEl ]. The nucleon case appears 
considerably more complicated. Both cases, however, can be analytically continued using 
the same method which yields ordinary one-dimensional integrals. 

Let us return to the expression for the neutral pion energy in a Euclidean electric field 
£. Setting the charge to unity e = 1, we have 
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The overall factor arises from the change to a dimensionless integration variable. Taking 
£ — > iE, where E is the electric field in Minkowski space, and defining x = m^/E, we have 



ioo ixs , . roo+ie 



iZ{-ix)=i ds— — — -1 = / 1. (C2) 

Jo Vsmhs / Jo+ie s vsms / 

In the last equality, we have added an ie prescription to render in the Minkowski space am- 
plitude well defined |l| . Because the integrand has only simple poles, we can evaluate using 
the Cauchy principle value for the real part and the sum of residues for the imaginary part 
(technically multiplied by one-half). This method is not general enough for our purposes, 
and we pursue a different route. 

By shifting the integration variable, all non-analyticities lie in the lower-half complex 
s-plane. We then consider the closed contour integral over the quarter circle in the first 
quadrant (positively oriented). This integral is zero by Cauchy; moreover, the integral 
around the quarter arc at infinity is zero by Jordan. Denoting the integrand by f(s), we 
have the equality 



poo rQ 

/ dsf(s + ie) + dsf(s) = 0. (C3) 

Jo J ioo 



Consequently we find 



oc 







= * ^ !,, (04) 



which gives a well-defined integral for the analytically continued function. Calculated nu- 
merically, the real part of Eq. (1C4I) agrees with the closed form expression given in Eq. (1531) . 
The imaginary part, 

$sm\iZ(-ix)] = - I" % cosixs) ( — ? , (C5) 

J s 2 \ smli s J 

we can compute in closed form with the aid of a Mittag-Leffler expansion 

4f-4-- l ) = 2 E^ z1 ^- (C6) 

s z Vsmhs / L — ' s z + n z ir z 

n=l 

Carrying out the proper time integration, leaves us with a sum 

%m[il(-ix)] = J^(-l) n - = - log (1 + e~ nx ) , (C7) 

n=l 



in agreement with the imaginary part of Eq. ( 1531) . 

The same analysis can be applied to analytically continue Eq. (11061) to determine the 
nucleon energy shift in an electric field. With x = m^/e\E\, the Minkowski space shift reads 

SE(A,-E>) = P~ XS < \l Eric [^^) [2 + scots] 

s 

+ - 



cossvcots 
2A s 3e sA2 / E 



'ttE cos s cot s 
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1 + 

E cot s 



(CI 



The integrand has poles and branch cuts. 10 Each of the singularities, by virtue of the ie 
prescription, is located in the lower-half plane. In the case of branch cuts, we must use the 
principal branch of each function to situate the cuts lie in the lower-half plane. We can again 
consider the integral over the quarter circle in the first quadrant to arrive at an expression 
for the analytic continuation. Defining y = xA 2 /m 2 , we have 



6E(A, —E 2 




2 iy 

e coth s Erfc 



e c ° th s Erfc 



Erfc 




[2 + s coth s] 



(C9) 



From this expression one can take the real and imaginary parts, although the results are 
quite cumbersome to display. Additionally useful are expressions for the error function of 
complex argument. For a and b real, we have 



9?e 



3m 



Erfc(a + ib) = Erfc(a) - 4=e~ a2 J^ 1 )™ 

* n=0 



Erfc(a + ib) 



n=0 
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(2n+ 1)! 



(2n + 2)\ 
H 2n (a,), 



H: 



2n+l( a J 



(CIO) 
(Cll) 



where H n (a) are Hermite polynomials. Simpler expressions result in the case when A = 0, 
and these have been quoted in the main text. 
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